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Abstract 

A measurement of the production cross-section for Z bosons that decay to muons 
is presented. The data were recorded by the LHCb detector during pp collisions 
at a centre-of-mass energy of 7 TeV, and correspond to an integrated luminosity 
of 1.0 fb _1 . The cross-section is measured for muons in the pseudorapidity range 
2.0 < r] < 4.5 with transverse momenta px > 20GeV/c. The dimuon mass is 
restricted to 60 < < 120GeV/c 2 . The measured cross-section is 

= ( 76 -° ± 0.3 A 0.5 A 1.0 ± 1.3) pb 

where the uncertainties are due to the sample size, systematic effects, the beam energy 
and the luminosity. This result is in good agreement with theoretical predictions 
at next-to-next-to-leading order in perturbative quantum chromodynamics. The 
cross-section is also measured differentially as a function of kinematic variables of the 
Z boson. Ratios of the production cross-sections of electroweak bosons are presented 
using updated LHCb measurements of W boson production. A precise test of the 
Standard Model is provided by the measurement of the ratio 

erw+^^+v^ + = 20 63 ± o.o9 ± 0.12 ± 0.05, 

VZ-s^+pi- 

where the uncertainty due to luminosity cancels. 
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1 Introduction 


Measurements of the total and differential cross-sections for the production of Z bosons in pp 
collisions test the Standard Model (SM) and provide constraints on parton density functions 
(PDFs) of the proton[] Theoretical predictions for these cross-sections are available at next- 
to-next-to-leading order (NNLO) in perturbative quantum chromodynamics (pQCD) [l]-j5]. 
The dominant uncertainty on these predictions reflects the uncertainties on the PDFs, 
which vary as functions of the kinematic variables studied. The forward acceptance of the 
LHCb detector allows the PDFs to be constrained at Bjorken-x values down to 10“ 4 |6|. 
Ratios of the W and Z cross-sections provide precise tests of the SM as the sensitivity 
to the PDFs in the theoretical calculations is reduced and many of the experimental 
uncertainties cancel. 

LHCb has measured the Z boson production cross-section at \/s= 7 TeV using decays 
to muon pairs in a data set corresponding to 37 pb _1 [7], and using electron [8] and tau 
lepton 191 pairs in a data set of 1.0 fb _1 . Production cross-sections of W bosons and 
the W + /W~ cross-section ratio have been measured in the muon channel 
1.0 fb” 1 data set. Similar measurements have also been performed by the ATLAS 
CMS 


' with the 
and 
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12 collaborations. 


The analysis described here is an update of the one described in Ref. |7|, using a total 
integrated luminosity of about 1.0 fb _1 . This increases statistical precision and allows 
better control of systematic uncertainties, with the result that the total uncertainties on 
the measurements are significantly reduced. Measurements are performed for muons with 
transverse momentum pt > 20 GeV/c and pseudorapidity in the range 2.0 < // < 4.5. In 
the case of Z boson measurements, the invariant mass of the two muons is required to be 
in the range 60 < M m +^- < 120 GeV/c 2 . These kinematic requirements define the fiducial 
region of the measurement and in this article are referred to as the fiducial requirements. 
Total cross-sections are presented as well as differential cross-sections as functions of the 
Z boson rapidity yz, Pt,z and (p* z . Here (f>* z is defined as 
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_ tan (0acop/2) 

z cosh(A?//2)' 


( 1 ) 


The angle </> acop = n — |A</>| depends on the difference A0 in azimuthal angle between the 
two muon momenta, while the difference between their pseudorapidities is denoted by Arj. 

The W boson cross-sections given in Ref. 1101 are re-evaluated using a more precise 
determination of the event trigger efficiency. The cross-sections are presented as a function 
of the 7] of the muon from the W boson decay. The values presented here supersede those 
of Ref. [101. 

This paper is organised as follows: Section [2] describes the LHCb detector; Sections [3] 
and [f] detail the selection of Z boson candidates, the Z boson cross-section definition and 
relevant sources of systematic uncertainty; Section [5] presents the results and Section [6] 
concludes the paper. Appendices [A] and [B] provide tables of differential cross-sections and 
correlations between these measurements. 

throughout this article Z represents both resonant production of Z bosons and off-mass-shell photons. 


1 








2 Detector and data set 


The LHCb detector [14|15] is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high-precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region |16jj, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of 
silicon-strip detectors and straw drift tubes 117 placed downstream of the magnet. The 
tracking system provides a measurement of momentum, p, of charged particles with a 
relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200GeV/c. The 
minimum distance of a track to a primary vertex, the impact parameter, is measured 
with a resolution of (15 + 29/px) b m , where px is the component of the momentum 
transverse to the beam, in GeV/c. Different types of charged hadrons are distinguished 
using information from two ring-imaging Cherenkov detectors 18 . Photons, electrons and 


hadrons are identified by a calorimeter system consisting of a scintillating-pad detector 
(SPD), preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. 
Muons are identified by a system composed of alternating layers of iron and multiwire 

The online event selection is performed by a trigger 


proportional chambers 19 


20 , which 


consists of a hardware stage, based on information from the calorimeter and muon systems, 
followed by a software stage, which applies a full event reconstruction. A requirement that 
prevents events with high occupancy from dominating the processing time of the software 
trigger is also applied. This is referred to as the global event cut (GEC) in this article. 

The measurement presented here is based on pp collision data collected at a centre- 
of-mass energy of 7 TeV. The integrated luminosity amounts to 975 ± 17 pb _1 . The 
absolute luminosity scale was measured during dedicated data taking periods, using both 
Van der Meer scans |21j and beam-gas imaging methods 22 . Both methods give similar 


results, which are combined to give the final luminosity estimate with an uncertainty of 
1.7% 


23 . This analysis uses the same data set as in Ref. 110 


Several samples of simulated data are produced to estimate contributions from back¬ 
ground processes, to cross-check efficiencies and to unfold data for detector-related effects. 

configured as in Ref. 


The PYTHIA generator 24,25 


26 with the CTEQ6L1 27,28 


parameterisation for the PDFs, is used to simulate 55, cc, WW, ti and Z production. All 
generated events are passed through a detector simulation based on Geant4 1291 , followed 
by LHCb-specific trigger emulation and event reconstruction. 

The results of the analysis are compared to theoretical predictions calculated with the 
30 31 generator at NNLO for the PDF sets ABM12 32 , CT10 33 , HERA1.5 
and NNPDF3.0 


Fewz 

JR09 


35 


MSTW08 36 


34 


37 . Comparisons are also made to the 


38-40 and Powheg 41 generators configured with the CT10 PDF set. ResBos 


ResBos 

includes an approximate NNLO calculation, plus a next-to-next-to-leading logarithm 
approximation for the resummation of the soft gluon radiation. Powheg provides a 
next-to-leading order (NLO) calculation interfaced to a parton shower, in this case per¬ 
formed by Herwig 142, 431. The results are also compared to the predictions from 
MC@NLO 44,45 , which is interfaced with different generators to simulate the parton 
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shower. Parton showering is performed using HERWIG 42,43 


the root mean-square-deviation of the intrinsic fey, and HERWIRI 46-48 

DGLAP-CS 


with different values for 
which is based 


on infrared-improved 49 DGLAP-CS 50-55 theory. All calculations are performed with 
the renormalisation and factorisation scales set to the electroweak boson mass. Scale 
uncertainties are estimated by varying these scales by factors of two around the boson 
mass 1561. Total uncertainties correspond to the PDF and a s uncertainties at 68.3% 
confidence level and scale uncertainties, added in quadrature. 


3 Event selection 

Events considered in this analysis are selected by the muon trigger. At the hardware stage, 
this trigger requires a muon with px > 1.5GeV/c and imposes an upper limit of 600 hits 
in the SPD sub-detector. At the software stage, a muon with px > lOGeV/c is required. 
The muon track must also satisfy additional track quality criteria. Candidate events are 
selected by requiring a pair of well-reconstructed particles of opposite charge, identified as 
muons, that also pass the fiducial requirements [7|. In total, 58466 Z boson candidates 
are selected and their invariant mass distribution is shown in Fig. [l] 



[GeV/c2] 


Figure 1: Invariant mass of dimuon candidates. 

The background contamination in the candidate sample is low. Five background sources 
are investigated: decays of heavy flavour hadrons, hadron misidentification, Z —y t + t~ 
decays, tt and W + W~ production. Unlike muons from signal, muons arising from decays 
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of heavy flavour hadrons are neither directly produced at the primary interaction vertex, 
nor are they isolated particles. Using the techniques from Ref. |7|, the contribution from 
this background is estimated from the data to be 227 ± 32 events, which amounts to 
0.4% of the candidate sample. The contribution from hadrons that decay in flight or 
have sufficient energy to traverse the calorimeters and be detected in the muon stations 
is studied in randomly triggered data, as described in Ref. [7], and determined to be 
116 ± 45 events, which is 0.2% of the candidate sample. Other electroweak and QCD 
backgrounds are estimated using PYTHIA simulation 25 and normalised to the measured 
total cross-sections for Z —y t + t~ decays 57 58 , tt 59,60 and W + W~ 61,62 production. 


The estimate for these sources is 66 ± 6 events, or 0.1% of the candidate sample. In total, 
the background is estimated to be 409 ± 56 events, or 0.7% of the candidate sample. 

The purity, defined to be the ratio of signal to total candidate events, is p=0.993±0.002. 
It is assumed to be constant as a function of y z , Pt,z and (j>* z . A systematic uncertainty, 
discussed later, is assigned to allow for possible inaccuracies in this assumption. 


4 Cross-section measurement 


Cross-sections are quoted in the kinematic range defined by the measurement and are 
corrected for quantum electro dynamic (QED) final-state radiation (FSR) in order to 
provide a consistent comparison with NLO and NNLO QCD predictions. No corrections 
are applied for initial-state radiation, electroweak effects, nor their interplay with QED 
effects. The cross-section in a given bin i of y z , Pt,z or 4>*zi with both final-state muons 
inside the fiducial region, is measured as 


& z —DJ 



£ £gec(* 


( 2 ) 


The indices i and j run over the bins of the variable under study. The index k runs over 
the candidates contributing to bin j. The total muon reconstruction efficiency for an event 
is given by e(r /£ , 77 ^ ), which is dependent on the pseudorapidity of the two muons. The 
matrix U corrects the data for bin migrations due to detector resolution effects. It is 


determined using an unfolding procedure, which is described in Section |4.4| The efficiency 
of the requirement on the number of SPD hits in the hardware trigger is denoted by 
£gec- The correction factors for QED final-state radiation are denoted by /fsr(*) and are 
determined for each bin. The integrated luminosity is denoted by £. Though not entering 
the expression for the cross-section, an uncertainty due to the beam energy is assigned to 
all cross-sections. More detail on these individual components is given below. Once the 
binned cross-sections are determined, they are summed to give the total cross-section 


Pz^/i+fi- 


'y ^ (*)■ 


(3) 


The most precise estimate of the total cross-section is obtained by summing the cross- 
sections as a function of rapidity, where uncertainties due to data unfolding are negligible. 
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4.1 Muon reconstruction efficiencies 


The data are corrected for efficiency losses due to track reconstruction, muon identification, 
and trigger requirements. All efficiencies are determined from data using the techniques 
detailed in Refs. |7,8], where the track reconstruction, muon identification, and muon 
trigger efficiencies are obtained using tag-and-probe methods on the Z resonance. The 
tag and probe tracks are required to satisfy the fiducial requirements. The tag must 
be identified as a muon and be consistent with triggering the event, while the probe is 
defined so that it is unbiased by the requirement for which the efficiency is being measured. 
The efficiency is studied as a function of several variables, which describe both the muon 
kinematics and the detector occupancy. In this analysis the efficiency as a function of 
muon r) is used. The efficiency in each bin of r) is defined as the fraction of tag-and-probe 
candidate events where the probe satisfies a track reconstruction, identification or trigger 
requirement. 

The tracking efficiency is determined using probe tracks that are reconstructed by 
combining hits from the muon stations and the large-area silicon-strip detector. The 
efficiency depends on r] and varies between 89.5% and 98.5% with uncertainties between 
0.4% and 1.9%. 

The muon identification efficiency is determined using probe tracks that are recon¬ 
structed without using the muon system. The efficiency depends on t] and varies between 
91.3% and 99.2% with uncertainties between 0.1% and 0.9%. 

The single-muon trigger efficiency is determined using reconstructed muons as probes. 
The efficiency depends on ij and varies between 71.6% and 82.0% with uncertainties 
between 0.5% and 1.2%. Since only one muon candidate is required for the event to pass 
the trigger requirements, the overall trigger efficiency for the analysis is about 95%. 

The efficiency to reconstruct any given event is taken to be the product of the three 
individual efficiencies and determined on an event-by-event basis as a function of muon 7 ], 


e(/i + ,/i ) = 


~ £ trk ' £ trk ' £ id 


'id 



+ £' 


trg 



(4) 


In Equation |4j the efficiency e is written explicitly in terms of the muon tracking (e tr k), 
identification (e id ) and trigger (£ trg ) efficiencies. The average reconstruction efficiency 
for the analysis is about 85%. Effects that correlate the efficiency of the two muons are 
considered, but these are negligible at the current level of precision. 


4.2 GEC efficiency 

The GEC efficiency is the efficiency of the SPD multiplicity limit at 600 hits in the muon 
trigger. This efficiency is evaluated from data using two independent methods. The 
first exploits the fact that the SPD multiplicities of single pp interactions are always 
below the 600 hit threshold. The expected SPD multiplicity distribution of signal events 
is constructed by adding the multiplicities of signal events in single pp interactions to 
the multiplicities of randomly triggered events, as in Ref. \7\. The convolution of the 
distributions extends to values above 600 hits, and the fraction of events that the trigger 
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rejects can be determined. The second method consists of fitting the SPD multiplicity 
distribution and extrapolating the fit function to determine the fraction of events that are 
rejected, as in Ref. [8i|. Both methods give consistent results and £gec = (94.0 ± 0.2)% 
is used in this analysis. The central value is the estimate from the first method, while 
the difference between the two estimates contributes to the uncertainty. This efficiency 
depends linearly on y z with about 2% variation across the full range. A weaker dependence 
on both the p? } z and cp* z is also observed. Corrections for these effects are made. 


4.3 Final-state radiation 

The FSR correction is taken to be the mean of the corrections calculated with 

. The corrections are tabulated in Appendix [A] and 

are on average about 2.5%. 


Herwig++ 1631 and Pythia8 25 


4.4 Unfolding detector response 

To correct for detector resolution effects, an unfolding is performed (matrix U of Equation [2]) 
using LHCb simulation and the RooUNFOLD [64] software package. The momentum 
resolution in the simulation is calibrated to the data. The data are then unfolded using the 
iterative Bayesian approach proposed in Ref. 65. Other unfolding techniques 66,67 give 


similar results. Additionally, all unfolding methods are tested for model dependence using 


underlying distributions from leading order PYTHIA 24 25 , leading order HERWIG+ + 
as well as NLO Powheg 


63 


|41|68,69] showered with both Pythia and Herwig using the 
Powheg matching scheme. The correction is on average about 2% as a function of Pt,z, 
while it is significantly less as a function of (j)* z . Only the px,z and (p* z distributions are 
unfolded. Since y z is well measured, no unfolding is performed and U is the identity 
matrix in this case. 


4.5 Systematic uncertainties 

Sources of systematic uncertainty and their effect on the total cross-section measurement 
are summarised in Table [T|p] The measured cross-sections as a function of px,z and cf z 
have additional systematic uncertainties due to unfolding. 

The systematic uncertainty associated with the trigger, identification and tracking 
efficiencies is determined by re-evaluating all cross-sections with the values of the individual 
efficiencies increased or decreased by one standard deviation. The full covariance matrix 
of the differential cross-section measurements is evaluated in this way for each source of 
uncertainty separately. The covariance matrices for each source are added and the diagonal 

2 Many of the systematic uncertainties quoted here have a statistical component. The statistical 
uncertainty quoted on the measurement is due to the number of observed Z candidates. In the case 
of unfolded measurements, the statistical uncertainty is provided by the covariance matrix returned by 
RooUnfold. 
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Table 1: Contributions to the relative uncertainty on the total Z boson cross-section. 


Source 

Uncertainty (%) 

Statistical 

0.39 

Trigger efficiency 

0.07 

Identification efficiency 

0.23 

Tracking efficiency 

0.53 

FSR 

0.11 

Purity 

0.22 

GEC efficiency 

0.26 

Systematic 

0.68 

Beam energy 

1.25 

Luminosity 

1.72 

Total 

2.27 


elements of the result determine the total systematic uncertainty due to reconstruction 
efficiencies. These vary between 0.5 and 2.0% on the differential cross-section measurements. 

The systematic uncertainty on the FSR correction is the quadratic sum of two com¬ 
ponents. The first is due to the statistical precision of the PYTHIA and HERWIG+ + 
estimates and the second is half of the difference between their central values. The latter 
dominates, with the uncertainties on the differential cross-sections varying between 0.3 
and 3%. 

The systematic uncertainty on the purity is determined from the number of candidate 
and background events. In addition, an uncertainty based on the assumption that the 
purity is the same for all variables and bins of the analysis is evaluated by comparing to 
cross-section measurements using a binned purity, rather than a global one. The total 
uncertainties on the differential cross-section measurements due to variations in purity are 
typically less than 1%. 

The GEC efficiency is determined in each bin of yz, Pt,z and <p* z . The systematic 
uncertainty is the sum in quadrature of a component due to the available sample size 
in each bin and a component due to the 0.2% uncertainty on the integrated number, 
as determined in Section |4.2l This varies between 0.4 and 4% across the differential 
measurements. 

The systematic uncertainty due to unfolding is estimated by the differences between 
the differential cross-sections using Bayesian and matrix inversion unfolding techniques. 
The typical size is 1.5%. 

The measurement is specified at centre-of-mass energy y/s = 7 TeV. The beam energy, 
and consequently the centre-of-mass energy, is known to 0.65% 70 . The sensitivity of the 


cross-section to the centre-of-mass energy is studied using DYNNLO 71 and a systematic 
uncertainty of 1.25% is assigned. 


7 










5 Results 


5.1 Z boson production cross-section 

The measured rapidity distribution as shown in Fig. [2] is compared to the prediction from 
Fewz 30, 31 with six different PDF sets. To compare the shapes of the differential 


41 


43 


cross-sections, measurements and predictions are normalised to the total fiducial cross 
section. The normalised differential cross-sections are shown in Figs. li and |4j The 
measurements are compared to the predictions from ResBos 38 40 1 and POWHEG 
where events are interfaced with a parton shower that is simulated using Herwig [42 
The pt,z and (J)* z distributions are well described by ResBos and POWHEG, with the 
central values overestimating the data slightly at low (J)* z and underestimating slightly at 
high (f>* z . Comparisons to MC@NLO + Herwiri and MC@NLO + Herwig are shown 
in Figs. [5] and [6] Here HERWIG is configured with the root mean-square-deviation of the 
intrinsic kx distribution set to OGeV/c in one instance and 2.2GeV/c in another. The 
predictions straddle the measurement at low px,z and (f>* z . The high px,z and <p* z tails are 
underestimated. 

The total inclusive cross-section for Z — y p + /j,~ production for muons with 
Pt > 20GeV/c in the pseudorapidity region 2.0 < r) < 4.5 and the dimuon invariant 
mass range 60 < < 120 GeV/c 2 is measured to be 


uz—tu+u,- = (76.0 ± 0.3 ± 0.5 ± 1.0 ± 1.3) pb, 


where the first uncertainty is statistical, the second systematic, the third is due to the 
beam energy and the fourth is due to the luminosity. The upper plot of Fig. [7] shows 
agreement between this measurement and NNLO predictions given by Fewz configured 
with various PDF sets. The measurement also agrees with the measurements of the Z 
boson production cross-section performed in the electron |8| and tau lepton |9| channels 
but with a significantly smaller uncertainty. All binned cross-sections are detailed in 
Tables [3j [4] and [5] of Appendix [Aj 


5.2 Ratios of electroweak boson production cross-sections 

The cross-section ratios are defined for muons with px > 20 GeV/c, 2.0 < p < 4.5 and, in 
the case of the Z boson cross-section, a dimuon invariant mass between 60 and 120 GeV/c 2 . 
The ratio of W boson to Z boson production is defined as 


Rwz — 


°~ W +- s-Ai+1^ + 

TZ-^+At- 


(5) 


The separate ratios of W + and W to Z boson production cross-sections are defined as 


Rw ± z 




( 6 ) 
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while the W boson cross-section ratio is defined as 

R W = . (7) 

Many sources of systematic uncertainty cancel or are reduced in the ratios. As the data 
sets are identical, the largest single source of uncertainty on the individual cross-sections, 
due to the luminosity determination, is removed. The trigger used to select both samples 
is identical and most of the uncertainty on the determination of the trigger efficiency 
cancels. In particular, the GEC is common to both the W and Z boson analyses and it is 
expected that the size of the efficiency correction is similar for W and Z events. Cross¬ 
checks in data and simulation support this assumption with a precision of approximately 
0.3%, which is included as a systematic uncertainty. The GEC efficiency was determined 
for the previous W boson measurement 1101 to be (95.9 ± 1.1)%, whereas an improved 
precision of (94.0 ± 0.2)% is obtained in the current analysis. Consequently, the W boson 
cross-section results are updated to benefit from the more precise value. These results 
are listed in Tables iEI and [8] of Appendix [Aj along with the muon charge ratios and 
asymmetries, and supersede those in Ref. [l0;. The uncertainties on the tracking and muon 
identification partially cancel in the ratios of W and Z bosons. The uncertainty on the 
W + {yV~) cross-section due to beam energy is 1.06(0.91)% and most of this uncertainty also 
cancels in the ratios. The uncertainties on the purities of the W and Z boson selections 
are uncorrelated. The FSR uncertainties are taken to be uncorrelated. The sources of 
uncertainty contributing to the determination of the ratios are summarised in Table [2j 

Table 2: Contributions to the relative uncertainty on the electroweak boson cross-section ratios. 


Source Uncertainty (%) 



Rwz 

Rw+z 

Rw-Z 

Rw 

Statistical 

0.45 

0.48 

0.50 

0.38 

Trigger efficiency 

0.15 

0.16 

0.13 

0.07 

Identification efficiency 

0.12 

0.12 

0.12 

0.03 

Tracking efficiency 

0.24 

0.23 

0.26 

0.08 

FSR 

0.16 

0.21 

0.17 

0.21 

Purity 

0.41 

0.49 

0.55 

0.62 

GEC efficiency 

0.27 

0.28 

0.29 

0.18 

Systematic 

0.60 

0.67 

0.72 

0.69 

Beam energy 

0.26 

0.19 

0.34 

0.15 

Total 

0.79 

0.85 

0.94 

0.80 


The dominant uncertainties on the ratios are due to the purity and the size of the 
sample. The correlation coefficients used in the uncertainty calculations are tabulated in 
Appendix |Bj 
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The updated W + boson cross-section is 


= (878.0 ± 2.1 ± 6.7 ± 9.3 ± 15.0) pb, 

where the uncertainties are due to the sample size, systematic effects, the beam energy 
and the luminosity determination. The updated W~ boson cross-section is 


a w~ 




-p = (689.5 ±2.0 ±5.3 ±6.3 ±11.8) pb. 


These measurements are in good agreement with the predictions of NNLO pQCD, as 
shown in Fig. [7J Using the Z boson cross-section from Section |5T , elect.roweak boson 
cross-section measurements and theoretical predictions, with different parameterisations 
of the PDFs, are compared in Fig. [8j with contours corresponding to the 68.3% confidence 

level. 

The W to Z boson cross-section ratio is measured as 


R wz = 20.63 ± 0.09 ± 0.12 ± 0.05, 

where the first uncertainty is statistical, the second is systematic and the third is due to 
the beam energy. The charged W to Z boson cross-section ratios are measured as 

R w+Z = 11.56 ± 0.06 ± 0.08 ± 0.02, 

R w - Z = 9.07 ±0.05 ±0.07 ±0.03, 

while the W boson cross-section ratio is measured as 


R w = 1.274 ± 0.005 ± 0.009 ± 0.002. 

These measurements, as well as their predictions, are displayed in Fig. [9j For Rwz and 
Rw+z , the data are well described by HERA1.5 and JR09, while the values from CT10, 
MSTW08, NNPDF3.0 and ABM12 are larger than those measured. All PDF sets show good 
agreement for R w -z- As previously reported [l0|, all PDF sets except ABM12 show good 
agreement for Rw- The Rwz and Rw ratios are measured with a fractional uncertainty of 
0.8%, which is similar both to the precision due to the PDFs on the individual theoretical 
predictions and to the spread between the predictions. Considering the spread in the 
different predictions, the experimental measurements are in good agreement with SM 
predictions and can be used to improve the determination of the PDFs. 


6 Conclusions 

A measurement of the forward Z boson production cross-section at y/s = 7 TeV is presented, 
where the Z bosons are reconstructed in the decay Z —> y + . The total cross-section 
in the fiducial range of the selection is in agreement with NNLO pQCD calculations. 
Normalised differential cross-sections as a function of yz, 4>* z and pr,z are compared to 
the predictions of various generators. The increased precision on the determination of 
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the event trigger efficiency motivates a re-evaluation of the recently measured W boson 
production cross-section. These are presented here and supersede the values given in 
Ref. 10 . Combining the Z boson cross-section with updated W boson cross-sections 


measured in a similar fiducial volume allows for precision measurements of elect.roweak 
boson cross-section ratios. In particular, the W to Z boson ratio is determined with a 
relative precision of 0.8%. The measured ratios are consistent with SM predictions but 
are sensitive to the particular choice of PDF. Consequently, these results are expected to 
provide significant constraints on PDFs. 


Acknowledgements 

We thank B.F.L. Ward for providing MC@NLO + Herwig and MC@NLO + Herwiri 
predictions. We express our gratitude to our colleagues in the CERN accelerator depart¬ 
ments for the excellent performance of the LHC. We thank the technical and administrative 
staff at the LHCb institutes. We acknowledge support from CERN and from the national 
agencies: CAPES, CNPq, FAPERJ and FINEP (Brazil); NSFC (China); CNRS/IN2P3 
(France); BMBF, DFG, HGF and MPG (Germany); INFN (Italy); FOM and NWO (The 
Netherlands); MNiSW and NCN (Poland); MEN/IFA (Romania); MinES and FANO 
(Russia); MinECo (Spain); SNSF and SER (Switzerland); NASU (Ukraine); STFC (United 
Kingdom); NSF (USA). The Tierl computing centres are supported by IN2P3 (France), 
KIT and BMBF (Germany), INFN (Italy), NWO and SURF (The Netherlands), PIC 
(Spain), GridPP (United Kingdom). We are indebted to the communities behind the 
multiple open source software packages on which we depend. We are also thankful for 
the computing resources and the access to software R&D tools provided by Yandex LLC 
(Russia). Individual groups or members have received support from EPLANET, Marie 
Sklodowska-Curie Actions and ERC (European Union), Conseil general de Haute-Savoie, 
Labex ENIGMASS and OCEVU, Region Auvergne (France), RFBR (Russia), XuntaGal 
and GENCAT (Spain), Royal Society and Royal Commission for the Exhibition of 1851 
(United Kingdom). 


11 





Figure 2: (a) Differential cross-section as a function of yz compared with the prediction of 
Fewz configured with various PDF sets. Different predictions are displaced horizontally for 
visibility, (b) Normalised differential cross-section as a function of yz compared to the predictions 
of ResBos and Powheg + Herwig. The shaded (yellow) bands indicate the statistical and 
total uncertainties on the measurements, which are symmetric about the central value. 
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Figure 3: Normalised differential cross-section as a function of pr,z on (a) logarithmic and (b) 
linear scales. The shaded (yellow) bands indicate the statistical and total uncertainties on the 
measurements, which are symmetric about the central value. The measurements are compared 
to the predictions of ResBos and Powheg + Hekwig. 
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Figure 4: Normalised differential cross-section as a function of (jf z on (a) logarithmic and (b) 
linear scales. The shaded (yellow) bands indicate the statistical and total uncertainties on the 
measurements, which are symmetric about the central value. The measurements are compared 
to the predictions of ResBos and Powheg + Hekwig. 
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Figure 5: Normalised differential cross-section as a function of pr,z on (a) logarithmic and (b) 
linear scales. The shaded (yellow) bands indicate the statistical and total uncertainties on the 
measurements, which are symmetric about the central value. The measurements are compared to 
MC@NLO + Herwig (HW) and MC@NLO + Herwiri (HERWIRI). Herwig is configured 
with two choices of the root mean-square-deviation of the intrinsic kx distribution, 0 and 2.2 
GeV/c. 
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Figure 6: Normalised differential cross-section as a function of 4>* z on (a) logarithmic and (b) 
linear scales. The shaded (yellow) bands indicate the statistical and total uncertainties on the 
measurements, which are symmetric about the central value. The measurements are compared to 
MC@NLO + Herwig (HW) and MC@NLO + Herwiri (HERWIRI). Herwig is configured 
with two choices of the root mean-square-deviation of the intrinsic kx distribution, 0 and 2.2 
GeV/c. 
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LHCb, is =7 TeV o MSTW08 □ ABM 12 > 20 GeV/c 

T 

Data stat v NNPDF30 o HERA 15 2.0 < rf < 4.5 

Data tot # CT10 a JR09 Z: 60 < < 120 GeV/c 2 




Figure 7: LHCb measurements of electroweak boson production cross-sections compared to 
NNLO pQCD as implemented by the Fewz generator using various PDF sets. The shaded 
(yellow) bands indicate the statistical and total uncertainties on the measurements, which are 
symmetric about the central value. 
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Figure 8: Two dimensional plots of electroweak boson cross-sections compared to NNLO predic¬ 
tions for various parameterisations of the PDFs. The outer, shaded (yellow) ellipse corresponds to 
the total uncertainty on the measurements. The inner, shaded (purple) ellipse excludes the beam 
energy and luminosity uncertainties. The uncertainty on the theoretical predictions corresponds 
to the PDF uncertainty only. All ellipses correspond to uncertainties at 68.3% confidence level. 
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Figure 9: Ratios of electroweak boson production R\y+z> R\v~Zi Rwz, Rw, compared to various 
theoretical predictions. The shaded (yellow) bands indicate the statistical and total uncertainties 
on the measurements, which are symmetric about the central value. 
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Appendices 

A Cross-sections 


Table 3: Inclusive differential cross-sections for Z boson production as a function of y z • Uncer¬ 
tainties are due to the sample size, systematic effects, the beam energy and the luminosity. No 
candidates are observed in the 4.250 - 4.500 bin. 


Vz 


a z [pb] 


/fsr 


2.000 

- 2.125 

0.969 

± 

0.039 

± 

0.032 

± 

0.012 

± 

0.017 

1.050±0.020 

2.125 

- 2.250 

2.840 

± 

0.063 

± 

0.050 

± 

0.036 

± 

0.049 

1.032±0.008 

2.250 

- 2.375 

4.428 

± 

0.077 

± 

0.078 

± 

0.055 

± 

0.076 

1.027±0.006 

2.375 

- 2.500 

5.823 

± 

0.088 

± 

0.060 

± 

0.073 

± 

0.100 

1.026±0.004 

2.500 

- 2.625 

6.877 

± 

0.095 

± 

0.068 

± 

0.086 

± 

0.118 

1.025±0.004 

2.625 

- 2.750 

7.669 

± 

0.100 

± 

0.069 

± 

0.096 

± 

0.132 

1.026±0.004 

2.750 

- 2.875 

8.306 

± 

0.104 

± 

0.070 

± 

0.104 

± 

0.143 

1.026±0.003 

2.875 

- 3.000 

8.241 

± 

0.103 

± 

0.066 

± 

0.103 

± 

0.142 

1.025±0.003 

3.000 

- 3.125 

7.783 

± 

0.099 

± 

0.059 

± 

0.097 

± 

0.134 

1.026±0.003 

3.125 

- 3.250 

7.094 

± 

0.096 

± 

0.058 

± 

0.089 

± 

0.122 

1.028±0.004 

3.250 

- 3.375 

5.894 

± 

0.087 

± 

0.049 

± 

0.074 

± 

0.101 

1.026±0.004 

3.375 

- 3.500 

4.160 

± 

0.073 

± 

0.041 

± 

0.052 

± 

0.072 

1.027±0.005 

3.500 

- 3.625 

2.896 

± 

0.061 

± 

0.030 

± 

0.036 

± 

0.050 

1.026±0.005 

3.625 

- 3.750 

1.741 

± 

0.047 

± 

0.023 

± 

0.022 

± 

0.030 

1.021±0.007 

3.750 

- 3.875 

0.825 

± 

0.032 

± 

0.014 

± 

0.010 

± 

0.014 

1.025±0.010 

3.875 

- 4.000 

0.321 

± 

0.020 

± 

0.008 

± 

0.004 

± 

0.006 

1.011±0.015 

4.000 

4.250 

- 4.250 

- 4.500 

0.115 

± 

0.013 

± 

0.006 

± 

0.001 

± 

0.002 

1.018±0.033 
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Table 4: Inclusive differential cross-sections for Z boson production as a function of pr,z- 
Uncertainties are due to the sample size, systematic effects, the beam energy and the luminosity. 


Tr,z[ GeV/c] 

oz [pb] 

/fsr 

0.0 - 

2.2 

6.454 ± 0.105± 0.129 ± 0.081 ± 0.111 

1.090T0.006 

2.2 - 

3.4 

6.520 ± 0.106± 0.150 ± 0.081 ± 0.112 

1.080T0.004 

3.4 - 

4.6 

6.209 ± 0.102± 0.221 ± 0.078 ± 0.107 

1.063T0.004 

4.6 - 

5.8 

5.868 ± 0.099± 0.208 ± 0.073 ± 0.101 

1.049T0.004 

5.8 - 

7.2 

5.749 ± 0.098± 0.154 ± 0.072 ± 0.099 

1.034T0.004 

7.2 - 

8.7 

5.607 ± 0.098± 0.083 ± 0.070 ± 0.096 

1.021T0.004 

8.7- 

10.5 

5.637 ± 0.098± 0.054 ± 0.070 ± 0.097 

1.002T0.004 

10.5 - 

12.8 

5.524 ± 0.096± 0.081 ± 0.069 ± 0.095 

0.996T0.004 

12.8 - 

15.4 

5.158 ± 0.092± 0.067 ± 0.064 ± 0.089 

0.984T0.005 

15.4 - 

19.0 

4.963 ± 0.087± 0.053 ± 0.062 ± 0.085 

0.978T0.005 

19.0 - 

24.5 

5.517 ± 0.088± 0.055 ± 0.069 ± 0.095 

0.985T0.004 

24.5 - 

34.0 

5.465 ± 0.085± 0.067 ± 0.068 ± 0.094 

1.013T0.004 

34.0 - 

63.0 

5.789 ± 0.085± 0.076 ± 0.072 ± 0.100 

1.038T0.004 

63.0 - 

270.0 

1.516 ± 0.043± 0.044 ± 0.019 ± 0.026 

1.060T0.007 


Table 5: Inclusive differential cross-sections for Z boson production as a function of 4>* z . Uncer¬ 
tainties are due to the sample size, systematic effects, the beam energy and the luminosity. 



oz [pb] 

/fsr 

0.00 - 0.01 

8.549 ± 0.099 ± 0.088 ± 0.107 ± 0.147 

1.034T0.004 

0.01 - 0.02 

7.805 ± 0.096 ± 0.106 ± 0.098 ± 0.134 

1.035±0.003 

0.02 - 0.03 

7.051 ± 0.091 ± 0.083 ± 0.088 ± 0.121 

1.034T0.004 

0.03 - 0.05 

11.362 ± 0.114 ± 0.108 ± 0.142 ± 0.195 

1.029T0.003 

0.05 - 0.07 

8.124 ± 0.097 ± 0.120 ± 0.102 ± 0.140 

1.026T0.003 

0.07-0.10 

8.436 ± 0.097 ± 0.074 ± 0.105 ± 0.145 

1.021±0.003 

0.10 - 0.15 

8.611 ± 0.098 ± 0.131 ± 0.108 ± 0.148 

1.020T0.003 

0.15 - 0.20 

4.819 ± 0.073 ± 0.092 ± 0.060 ± 0.083 

1.018T0.004 

0.20 - 0.30 

5.206 ± 0.076 ± 0.058 ± 0.065 ± 0.090 

1.019T0.004 

0.30 - 0.40 

2.541 ± 0.054 ± 0.051 ± 0.032 ± 0.044 

1.022T0.006 

0.40 - 0.60 

2.018 ± 0.048 ± 0.060 ± 0.025 ± 0.035 

1.024T0.007 

0.60 - 0.80 

0.755 ± 0.029 ± 0.035 ± 0.009 ± 0.013 

1.029T0.011 

0.80 - 1.20 

0.457 ± 0.023 ± 0.018 ± 0.006 ± 0.008 

1.025T0.014 

1.20 - 2.00 

0.166 ± 0.014 ± 0.011 ± 0.002 ± 0.003 

1.030T0.023 

2.00 - 4.00 

0.045 ± 0.008 ± 0.017 ± 0.001 ± 0.001 

1.031±0.041 
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Table 7: W + to W~ boson production cross-section ratios as a function of muon r/. Uncertainties 
are due to the sample size, systematic effects and the beam energy. These supersede the results 
in Ref. 


10 


rf 

Rw 

2.00 - 

2.25 

1.730 

± 

0.018 

± 

0.030 

± 

0.003 

2.25 - 

2.50 

1.706 

± 

0.015 

± 

0.040 

± 

0.003 

2.50 - 

2.75 

1.606 

± 

0.014 

± 

0.021 

± 

0.002 

2.75 - 

3.00 

1.388 

± 

0.013 

± 

0.024 

± 

0.002 

3.00 - 

3.25 

1.169 

± 

0.012 

± 

0.021 

± 

0.002 

3.25 - 

3.50 

0.898 

± 

0.010 

± 

0.025 

± 

0.001 

3.50 - 

4.00 

0.626 

± 

0.007 

± 

0.006 

± 

0.001 

4.00 - 

4.50 

0.328 

± 

0.011 

± 

0.011 

± 

0.000 


Table 8: Lepton charge asymmetries as a function of muon ?/. Uncertainties are due to the sample 
size, systematic effects and the beam energy. These supersede the results in Ref. 110 . 





4* 

[% 

] 



2.00 - 

2.25 

26.74 

± 

0.48 

± 

0.82 

± 

0.07 

2.25 - 

2.50 

26.08 

± 

0.41 

± 

1.09 

± 

0.07 

2.50 - 

2.75 

23.25 

± 

0.42 

± 

0.60 

± 

0.07 

2.75 - 

3.00 

16.26 

± 

0.46 

± 

0.84 

± 

0.07 

3.00 - 

3.25 

7.81 

± 

0.50 

± 

0.90 

± 

0.08 

3.25 - 

3.50 

-5.37 

± 

0.57 

± 

1.35 

± 

0.08 

3.50 - 

4.00 

-23.04 

± 

0.52 

± 

0.49 

± 

0.07 

4.00 - 

4.50 

-50.65 

± 

1.22 

± 

1.30 

± 

0.06 
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B Correlation Matrices 
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